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ABSTRACT 

To derive a new Ha luminosity function and to understand the clustering properties of star-forming 
galaxies at z k, 0.24, we have made a narrow-band imaging survey for Ha emitting galaxies in the 
HST COSMOS 2 square degree field. We used the narrow-band filter NB816 (A c = 8150 A, AA = 120 
A) and sampled Ha emitters with £W b s (Ha + [Nil]) > 12 A in a redshift range between z = 0.233 
and z — 0.251 corresponding to a depth of 70 Mpc. We obtained 980 Ha emitting galaxies in a sky 
area of 5540 arcmin 2 , corresponding to a survey volume of 3.1 x 10 4 Mpc 3 . We derive a Ha luminosity 
function with a best-fit Schechter function parameter set of a = — 1.35~^q'\\, log</>* = — 2.65^Q'3g, 

and log£*(erg s" 1 ) = 41.94±g'||. The Ha luminosity density is 2.7t°ol x 1q39 er S s s_1 Mpc" 3 . After 
subtracting the AGN contribution (15 %) to the Ha luminosity density, the star formation rate density 
is evaluated as 1.8~^p'4 x 10~ 2 M Q yr _1 Mpc -3 . The angular two-point correlation function of Ha 
emitting galaxies of logL(Ha) > 39.8 is well fit by a power law form of w{6) = 0.013±g ;™ 2 6»- a88±0 - 03 , 
corresponding to the correlation function of £(r) = (r/1.9Mpc) -188 . We also find that the Ha emitters 
with higher Ha luminosity are more strongly clustered than those with lower luminosity. 
Subject headings: galaxies: distances and redshifts — galaxies: evolution — galaxies: luminosity 
function, mass function 
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1. INTRODUCTION 

It is important to understand when and where intense 
star formation occurred during the course of galaxy evo- 
lution. Although the star formation history in individual 
galaxies is interesting, a general trend of star formation 
in galaxies as a function of time (or redshift) also provides 
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important insights on the global star formation history 
as well as on the metal enrichment history in the uni- 
verse. Therefore, the star formation rate density (SFRD) 
is one of the important observables for our understand- 
ing of galaxy formation and evolution. In the last decade, 
many works have followed the pioneer work of Madau et 
al. (1996) which compiled the evolution of SFRD, psfr, 
as a function of redshift for the first time. The evolution 
of Psfr is now widely accepted as follows: psfr steeply 
increases from z ~ to z ~ 1, and seems to be con- 
stant between z ~ 2 and z ~ 5 and may decline beyond 
z ~ 5 (Hopkins 2004 and references therein; Giavalisco et 
al. 2004; Taniguchi et al. 2005; Bouwens & Illingworth 
2006). 

Recent observations by the Galaxy Evolution Explorer 
(GALEX) and the Spitzer Space Telescope have con- 
firmed that psfr increases from z ~ to z ~ 1 (e.g., 
Schiminovich et al. 2005; Le Floc'h et al. 2005). How- 
ever, their observations show that the IR luminosity den- 
sity evolves as (1 + z) 4 while the UV luminosity density 
evolves as (1 + z) 2 ' 5 . This may imply that extinction by 
dust and reradiation from dust becomes to play a more 
important role at higher redshift. One of the remaining 
problems in this field is a relation between star-formation 
activity and large-scale structure formation. To study 
this issue, wide-field deep surveys are important. 

There are several star formation rate (SFR) estima- 
tors, e.g., UV continuum, Ha emission, [On] emission, 
far- infrared (FIR) emission (Kennicutt 1998), and radio 
continuum (Condon 1992). Each estimator has both ad- 
vantage and disadvantage to estimate SFR. UV contin- 
uum and nebular emission lines are considered to be di- 
rect tracers of hot massive young stars. However, they 
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are often affected by dust obscuration. On the other 
hand, FIR and radio continuum are insensitive to dust 
obscuration. FIR emission is due to the dust heated by 
the general interstellar radiation field. If most of the 
bolometric luminosity of a galaxy absorbed by dust is 
radiated from young stars, as in the case of dusty star- 
bursts, the FIR luminosity is a good SFR estimator. For 
early- type galaxies, much of the FIR emission is consid- 
ered to be related to the old stars and the FIR emission 
is not a good SFR estimator (Sauvage & Thuan 1992; 
Kennicutt 1998). For star- forming galaxies, there is a 
tight radio-FIR correlation (Condon 1992). This rela- 
tion suggests that the radio continuum also provides a 
good SFR estimator. The radio continuum is considered 
to be dominated by synchrotron radiation from relativis- 
tic electrons which are accelerated in supernova remnants 
(SNRs) (Lequeux 1971; Kennicutt 1983a; Gavazzi, Coc- 
ito, & Vcttolani 1986). We note that the radio continuum 
emission of some galaxies is dominated by the AGN com- 
ponent, although such galaxies are distinguished from 
star-forming galaxies by using the tight radio-FIR cor- 
relation (Sopp & Alexander 1991; Condon 1992). The 
nearly linear radio-FIR correlation also suggests that ra- 
dio continuum is affected by the efficiency of cosmic-ray 
confinement, since the degree of dust attenuation be- 
comes larger for more luminous galaxies (Bell 2003). Al- 
though SFRs evaluated from different SFR estimators 
are consistent with each other within a factor of 3 if the 
appropriate correction is applied for each case (e.g., Hop- 
kins et al. 2003; Chariot & Longhetti 2001; Chariot ct 
al. 2002), samples selected with a different method may 
have different biases. For example, samples selected by 
an objective-prism imaging survey are biased toward the 
system with large equivalent width (e.g., Gallego et al. 
1995), while those selected by UV radiation are biased 
against heavily dusty galaxies (Meurer et al. 2006). To 
evaluate the true SFRD, it is important to correct the ob- 
tained SFR appropriately and to know probable biases 
for the sample selection. 

In this work, we use the Ha luminosity as a SFR es- 
timator. The Ha luminosity is directly connected to 
the ionizing photon production rate. There are two ap- 
proaches to measure Ha luminosities of galaxies. One 
is a spectroscopic survey and the other is a narrow-band 
imaging survey. Although spectroscopic observations tell 
us details of emission line properties, e.g., Balmer decre- 
ment, metallicity, and so on, it is difficult to obtain spec- 
tra of a large sample of faint galaxies. On the other 
hand, narrow-band imaging observations make it possible 
to measure an emission-line flux of galaxies over a wide 
field of view. Another advantage of narrow-band imag- 
ing is that aperture corrections dose not need to eval- 
uate the total flux of Ha emission. However, there are 
some shortcomings in this method: e.g., narrow-band fil- 
ter cannot separate Ha emission from [Nil] AA6548, 6583 
emission and we cannot evaluate the obscuration degree 
for each galaxy. Therefore, we must correct these effects 
statistically. Since the redshift coverage of emission-line 
galaxies discovered by the narrow-band imaging method 
is restricted, the survey volume of emission-line galaxies 
is small. Therefore, it is difficult to obtain a large sample 
of Ha emitters. If this is the case, brighter (i.e., rarer) 
Ha emitters could be missed in such an imaging survey. 
In order to study the Ha luminosity function unambigu- 



ously, we need a large sample of Ha emitters covering a 
wide range of Ha luminosity. On the other hand, this 
restriction allows us to investigate large-scale structures 
of emission-line galaxies (mostly, star-forming galaxies) 
at a concerned redshift slice. 

Motivated by this in part, we have carried out a 
narrow-band imaging survey of the HST COSMOS 
field centered at a(J2000)= 10 h 00 m 28.6 s and (5(J2000)= 
+02°12'21.0"; the Cosmic Evolution Survey (Scoville et 
al. 2007). Since this field covers 2 square degree, it is 
suitable for our purpose. Our optical narrow-band imag- 
ing observations of the HST COSMOS field have been 
made with the Suprime-Cam (Miyazaki et al. 2002) on 
the Subaru Telescope (Kaifu et al. 2000; lye et al. 2004). 
Since the Suprime-Cam consists of ten 2kx4k CCD chips 
and provides a very wide field of view (34' x 27'), this is 
suitable for any wide- field optical imaging surveys. In our 
observations, we used the narrow-passband filter, NB816, 
centered at 8150 A with the passband of AA = 120 A. 
Our NB816 imaging data are also used to search both for 
Lya emitters at z m 5.7 (Murayama et al. 2007) and for 
[Oil] emitters at z « 1.2 (Takahashi ct al. 2007). In this 
paper, we present our results on Ha emitters at z w 0.24 
in the HST COSMOS field. 

Throughout this paper, magnitudes are given in the 
AB system. We adopt a flat universe with f2 ma tter = 0.3, 
Q A = 0.7, and H = 70 km s" 1 Mpc" 1 . 

2. PHOTOMETRIC CATALOG 

In this analysis, we use the COSMOS official photo- 
metric redshift catalog which includes objects whose to- 
tal i magnitudes [i! or i*) arc brighter than 25. The 
catalog presents 3" diameter aperture magnitude of 
Subaru/Suprime-Cam B, V, r' , i', z', and NB816 18 . 
Details of the Suprime-Cam observations are given in 
Taniguchi et al. (2007). Details of the COSMOS official 
photometric redshift catalog is also described in Capak 
et al. (2007) and Mobasher et al. (2007). Since the 
accuracy of standard star calibration (±0.05 magnitude) 
is too large to obtain an accurate photometric redshift, 
Capak et al. (2007) re-calibrated the photometric zero- 
points for photometric redshift using the SEDs of galax- 
ies with spectroscopic redshift. Following the recommen- 
dation of Capak et al. (2007), we apply the zero-point 
correction to the photometric data in the official cata- 
log. The offset values are 0.189, 0.04, -0.040, -0.020, 
0.054, and -0.072 for B, V, r 1 , i', z', and AB816, re- 
spectively. The zero-point corrected limiting magnitudes 
are B = 27 A, V = 26.5, r> = 26.6, i' = 26.1, z' = 25.4, 
and NB816 = 25.6 for a 3a detection on a 3" diameter 
aperture. The catalog also includes 3" diameter aperture 
magnitude of CFHT i*. We use the CFHT i* magni- 
tude for bright galaxies with i' < 21 because such bright 
galaxies appear to be slightly affected by the saturation 
effect in i' obtained with Suprime-Cam. We also apply 
the Galactic extinction correction adopting the median 
value E(B - V) = 0.0195 (Capak et al. 2007) for all 
objects. A photometric correction for each band is as 
follows (see Table 8 of Capak et al. 2007): A B = 0.079, 

18 Our SDSS broad- band niters are designated as g+, r+, i+, 
and z+ in Capak et al. (2007) to distinguish from the original 
SDSS filters. Also, our B and V filters arc designated as Bj and 
Vj in Capak ct al. (2007) where J means Johnson and Cousins 
filter system used in Landolt (1992). 
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A v = 0.061, A r , = 0.050, A*/ = 
A N saw = 0.034, and Ai- = 0.037. 



3. RESULTS 



3.1. 



Selection of NB816-Excess Objects 

We select Ha emitter candidates using 3" diameter 
aperture magnitude in the official catalog. In order to 
select ./VB816-excess objects efficiently, we need magni- 
tude of frequency-matched continuum. Since the effec- 
tive frequency of the NB816 filter (367.8 THz) is different 
either from those of i' (394.9 THz) and z' (333.6 THz) 
filters, we newly make a frequency-matched continuum, 
"iz continuum" , using the following linear combination 
; fiz — 0.57 fi> + 0.43/ z / where fi> and f z * are the i' and 
z' flux densities, respectively. Its 3 er limiting magnitude 
is iz ~ 26.03 in a 3" diameter aperture. For the bright 
galaxies with i' < 21, "iz continuum" is calculated as 
fi Z = 0.57 fi* + 0.43/ z ', where fi* is the i* flux density, 
since i' magnitude is incorrect because of the saturation 
effect. 

Since we use the ACS catalog prepared for studying 
weak lensing (Leauthaud et al. 2007) to separate galax- 
ies from stars, our survey area is restricted to the area 
mapped in Igu band with Advanced Camera for Sur- 
veys (ACS) on HST. After subtracting the masked out 
area, the effective survey area is 5540 arcmin 2 . Since 
the covered redshift range is between 0.233 and 0.251 
(Az = 0.018) and the corresponding survey depth is 70 
Mpc, our effective survey volume is 3.1 x 10 4 Mpc 3 . 

We selected iVB816-excess objects using the following 
criteria: 

iz - NB816 > 0.1, 



and 



iz 



NB816 > 3a(iz - NB816), 



where 

3a(iz-NB816) 



-2.51og(W(/ 3CT1VB8I6 ) 2 + (/ 3CTi J 2 //. 

(3) 

In the calculation of 7>a(iz — NB816), we applied the 
Galactic extinction correction to the limiting magnitudes 
of i'- and z'-band. The former criterion corresponds 
EW bs > 12 A. This criterion is exactly same as that 
of Fujita et al. (2003) and similar to that of Tresse & 
Maddox (1998) [EW(Ha + [NII]) rost > 10 A]. Taking ac- 
count of the scatter of iz — NB81Q color, we added the 
latter criterion. These two criteria are shown by the solid 
and dashed lines, respectively, in FigureQ] As we will de- 
scribe in the next section, we use the broad-band colors of 
galaxies to separate Ha emitters from other emission-line 
galaxies. To avoid the ambiguity of broad-band colors, 
we select galaxies detected above 3ct in all bands. Finally, 
we find 6176 galaxies that satisfy the above criteria. 



0.037, A z i = 0.028, their broad-band color properties comparing observed 

colors of our 6176 emitters with model ones that are es- 
timated by using the model spectral energy distribution 
derived by Coleman, Wu, & Weedman (1980). In Fig- 
ures [2] & [3l we show the B — V vs. V — r' and B — r' 
vs. i' — z' color-color diagram of the 6176 sources and 
the loci of model galaxies. Then we find that Ha emit- 
ters at z « 0.24 can be selected by adopting the fol- 
lowing three criteria; (1) B — V > 2(V - r') ~ 0.2, (2) 
B-r' > 5(z'-z')-1.3, and (3) B-r' > 0.7(i'-z') + 0.4. 
We can clearly distinguish Ha emitters from [Oin] or H/3 
emitters using the first criterion. We can also distin- 
guish Ha emitters from [On] emitters using the second 
and third criteria. We have checked the validity of our 
photometric selection criteria using both the photomet- 
ric data and spectroscopic redshifts of galaxies in the 
GOODS-N region (Cowie et al. 2004). Galaxies with 
redshifts corresponding to our Ha, [Oin], H/3, and [On] 
emitters are separately plotted in Figs. 2 & 3. It is 
shown that our criteria can separate well Ha emitters 
from [Oin], H/3, and [On] emitters. These criteria give 
us a sample of 981 Ha emitting galaxy candidates. The 
properties of GOODS-N galaxies presented in Figs. 2 
and 3 suggest that there is few contamination in our Ha 
emitter sample. 

3.3. Ha Luminosity 

As we mentioned in section 1, one of the advantages of 
narrow-band imaging is to measure the total flux of Ha 
emission directly without any aperture correction. To 
derive the total Ha flux, we have used the total flux of i' 

(1) (or i*), z', and NB816 using public images. 19 Our pro- 
cedure is the same as that given in Capak et al. (2007); 
MAG_AUTO in SExtractor (Bertin & Arnouts 1996). 

(2) Because of the contamination of the foreground galax- 
ies, one galaxy has a negative value of iz — NB816 based 

Sn the total magnitudes. We do not use this object in 
Irther analysis. Therefore, our final sample contains 980 



Ha emitters. 

Adopting the same method as that used by Pascual et 
al. (2001), we express the flux density in each filter band 
as the sum of the line flux, F^, and the continuum flux 
density, / c : 



/nb — fc 



fv =fc + -r^ 



ANB ' 
Fl_ 
Ai' 



(4) 
(5) 



and 

where ANB 
NB816 and i 
expressed as 



/*' - fc, (6) 

and Ai' are the effective bandwidths of 
, respectively. The iz continuum, /. 



is 



3.2. 



0.24 



Selection of NB816-Excess Objects at z 

The emission-line galaxy candidates selected above in- 
clude not only Ha emitters at z = 0.24 but also possibly 
[Oin] emitters at z = 0.63, or H/3 emitters at z = 0.68, 
or [On] emitters at z = 1.19 (Tresse et al. 1999; Ken- 
nicutt 1992b). We also note here that the narrowband 
filter passband is too wide to separate [Nn]AA6548, 6583 
from Ha. 

In order to distinguish Ha emitters at z « 0.24 from 
emission-line objects at other redshifts, we investigate 



fi. = 0.57 U + 0.43/* =fc + 0.57-^ 

Ai' 



(7) 



Using equation 0] and [71 the line flux Fl is calculated by 

/NB — fiz 



Ft, = ANB- 



(8) 



1 - 0.57{ANB/Ai') 

The line flux evaluated above includes both Ha and 
[Nil] AA6548, 6583 emission since the narrow-band filter 

19 http://irsa.ipac.caltech.edu/data/COSMOS/ 



4 



Shioya et al. 



cannot separate the contribution of these lines. The flux 
of Ha emission line is also affected by the internal extinc- 
tion. Therefore, we have to correct the contamination 
of [Nn]AA6548, 6583 emission and the internal extinction 
An a . Although several correction methods have been 
proposed (e.g., Kennicutt 1992a; Gallego et al. 1997; 
Tresse et al. 1994; Helmboldt et al. 2004 for [Nil] con- 
tamination: Kennicutt 1983b; Niklas et al. 1997; Ken- 
nicutt 1998; Hopkins et al. 2001; Afonso et al. 2003 
for Ana), there is few study which gives both corrections 
based on a single sample of galaxies. Helmboldt et al. 
(2004) have derived the relation between [Nil] /Ha and 
Mr and that between A^ a and Mr based on the data 
of the Nearby Field Galaxy Survey (Jansen et al. 2000a, 
2000b). We therefore adopt their relations to correct the 
[Nil] contamination and Aa a . After correcting to the 
AB magnitude system (Meurer et al. 2006), the relation 
between [Nil] /Ha and Mr is 

log ™ 6 583 = -0.13M R - 3.30, (9) 



with 



where 



^6583 



F 

_ [NII]6583A 



(10) 



and that between Ah q and Mr is 

log A Ha = -0.124^-2.47. (11) 

To derive Mr used in equations © & (fTTjl for each 
galaxy, we assume that the redshift of the galaxy is 
z = 0.242. We have also calculated fc-correction using the 
average SED of Coleman et al. (1980)'s Sbc and Irr. Tak- 
ing account of the luminosity distance and /c-correction 
(average value of Scd and Irr), Mr is calculated from 
r'-band total magnitude, r', as Mr = r' — 40.90. 

In addition to the above corrections, we also apply a 
statistical correction (21%; the average value of flux de- 
crease due to the filter transmission) to the measured flux 
because the filter transmission function is not square in 
shape (Fujita et al. 2003). Note that this value is slightly 
different from that (28 %) used in Fujita et al. (2003). 
Our new value is re-estimated by using the latest filter 
response function. The Ha flux is given by: 



x 10 4Ah ° x 1.21. 



(12) 

Finally the Ha luminosity is estimated by L(Ha) = 
4:Trd^F cor (Ra). In this procedure, we assume that all 
the Ha emitters are located at z = 0.242 that is the 
redshift corresponding to the central wavelength of our 
NB816 filter. Therefore, the luminosity distance is set to 
be d L = 1213 Mpc. 

We summarize the total magnitude of i' , z', NB816, 
and iz and the color excess of iz — NB816 for our Ha 
emission-line galaxy candidates in Table 1. Table 1 also 
includes logFL, logF cor (Ha), and logL(Ha). 

4. DISCUSSION 

4.1. Luminosity function of Ha emitters 

Figure 0] shows the Ha luminosity function (LF) at 
z pa 0.24 for our Ha emitter sample. The Ha LF is 
constructed by the relation 

^o g L t ) = -l-y^ (13) 



AlogL 



logij -logLi) < - AlogL, 



(14) 



V, 



where Vj is the volume of the narrow band slice in the 
range of redshift covered by the filter. We have used 
AlogL(Ha) = 0.2. If the shape of the filter response 
is square, our survey volume is 3.1 x 10 4 Mpc 3 . How- 
ever, effective survey volume is affected by the shape of 
filter transmission curve. For example, since the trans- 
mission at 8092 A is a half of the peak value, the color 
excess, iz — iV£>816, of an Ha emitter at z = 0.233 with 
EW(Ra + [Nil]) = 12A is observed as 0.05 which does 
not satisfy our selection criterion, iz — NB816 > 0.1. 
Taking account of the filter shape in the computation of 
the volume, the correction can be as large as 23 % for 
the faintest galaxies. Adopting the Schechter function 
form (Schechter 1976), we obtain the following best-fit 
parameters for our Ha emitters with L(Ha) > 10 39 8 
ergs s" 1 ; a = -1.35±g;^, log<^ = -2.65±g;§|, and 
logL*(erg s- 1 ) = 41.94l°j| (black solid line). 

Together with our Ha LF, Figure [4] shows Ha LFs of 
previous studies in which Ha emitters at z < 0.3 are 
investigated; Tresse & Maddox (1998) [which is char- 
acterized by a = -1.35, 0* = 10~ 256 Mpc" 3 , and 
L* = 10 41 ' 92 ergs s ; note that these parameters were 
converted by Hopkins (2004) to those of our adopted cos- 
mology], Fujita et al. (2003), Hippelein et al. (2003) and 
Ly et al. (2007). Fujita et al. (2003), Ly et al. (2007), 
and this work are based on the NB816 imaging obtained 
with the Subaru Telescope. Tresse & Maddox (1998) 
is based on the Canada-France Redshift Survey (CFRS) 
and Hippelein et al. (2003) is based on the Calar Alto 
Deep Imaging Survey (CADIS). 

First, we compare our Ha LF with that derived by Ly 
et al. (2007). Their best-fit Schechter function parame- 
ters (a = —1.71, log0* = —3.7, logL* = 42.2) are quite 
different from those of our Ha LF. However we note that 
the data points between logL(Ha) ~ 39.5 and ~ 41.0, 
shown in Fig. 10b of Ly et al. (2007), are quite simi- 
lar to our results (Figure |4]) . We therefore consider that 
the Ha LF of Ly et al. (2007) itself is basically con- 
sistent with ours except the brightest point. The differ- 
ence of Schechter parameters between ours and Ly et al's 
may arise from the data points of the brightest and the 
faintest ones, especially the brightest one. Since the field 
of view of the COSMOS is about an order wider than 
that of the SDF, we consider that our Ha LF is more 
accurate by determined than that of Ly et al. (2007) at 
the bright end. 

Second, we compare our Ha LF with the other Ha 
LFs. Although our Ha LF is similar to those of Tresse & 
Maddox (1998) and Hippelein et al. (2003), the Ha LF of 
Fujita et al. (2003) shows a steeper faint-end slope and 
a higher number density for the same luminosity than 
ours. These differences may be attributed to the follow- 
ing different source selection procedures: (1) Fujita et 
al. (2003) used their iVB<?l(}-selected galaxies while we 
used j'-selected galaxies, Tresse & Maddox (1998) used I- 
selected Canada-France Redshift Survey (CFRS) galax- 
ies, and Hippelein et al. (2003) used Fabry-Perot images 
for pre-selection of emission-line galaxies. As Fujita et al. 
(2003) demonstrated, samples based on a broad-band se- 
lected catalog are biased against galaxies with faint con- 
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tinuum. (2) Fujita et al. (2003) used their B — Rc vs. 
Rc — Ic color - color diagram to isolate Ha emitters from 
other low-z emitters at different redshifts. However, we 
find that there are possible contaminations of [Oin] emit- 
ters if one uses the B — Rc vs. Rc — Ic diagram, because 
of the small difference between Ha and [Oni] emitters on 
that color - color diagram. On the other hand, we used 
B — V vs. V — r' to isolate Ha emitters from [Oiii] emit- 
ters. Due to the large separation between Ha emitters 
and [Oiii] emitters on the B — V vs. V — r' diagram, 
we can reduce the contamination of [Oiii] emitters. (3) 
Fujita et al. (2003) used population synthesis model GIS- 
SEL96 (Bruzual & Chariot 1993) to determine the cri- 
teria for selecting Ha emitters. To check the validity 
of the criterion, we compare colors of GOODS-N galax- 
ies at z ~ 0.24, 0.63, 0.68 & 1.19 with model colors at 
corresponding redshifts based on GISSEL96 (Figure [5]). 
Unfortunately, the predicted colors are slightly different 
from those of observed galaxies. We therefore redeter- 
mined the selection criteria using the SED of Coleman, 
Wu, & Weedman (1980) as 

(B - R c ) > 2.5(i? c - Ic) + 0-2. 

If we adopt this revised criterion, the number of Ha emit- 
ters in the Fujita et al. (2003) is reduced by about 20 % 
(Figure[5]). This is one reason why the number density of 
Ha emitters in Fujita et al. (2003) is higher than other 
surveys. Recently, Ly et al. (2007) pointed out that the 
fraction of [Oin] emitters in the Ha emitter sample of 
Fujita et al. (2003) may be about 50 % using the Hawaii 
HDF-N sources with redshifts observed as iVi3816-excess 
objects. The Ha LF of Fujita et al. (2003) reduced by 
50 % appears to be quite similar to our Ha LF. 

4.2. Luminosity density and star formation rate density 
By integrating the luminosity function, i.e., 

C(Ra)= $(L)LdL = T(a + 2)faL„ (15) 
Jo 

we obtain a total Ha luminosity density of 2.7lgg x 
ergs s _1 Mpc -3 at z « 0.24 from our best fit LF. The star 
formation rate is estimated from the Ha luminosity using 
the relation SFR = 7.9 x 10~ 42 L(Ha) Mgyr" 1 , where 
L(Ha) is in units of ergs s _1 (Kennicutt 1998). Using 
this relation, the Ha luminosity density can be translated 
into the SFR density of p SFR ~ 2.1±^ x 10 _2 M Q yr" 1 
Mpc -3 . 

However, not all the Ha luminosity is produced by star 
formation, because active galactic nuclei (AGNs) can also 
contribute to the Ha luminosity. For example, previ- 
ous studies obtained the following estimates; 8-17% of 
the galaxies in the CFRS low-z sample (Tresse et al. 
1996), 8% in the Universidad Complutense de Madrid 
(UCM) survey of local Ha emission line galaxies (Gal- 
lego et al. 1995), and 17-28% in the 15R survey (Carter 
et al. 2001). Recently, Hao et al. (2005) obtained an 
Ha luminosity function of active galactic nuclei based on 
the sample of the Sloan Digital Sky Survey within a red- 
shift range of < z < 0.15. The Ha luminosity density 
calculated from Schechter function parameters which are 
shown in the paper is 1.1 x 10 38 erg s _1 Mpc -3 (with no 
reddening correction). Taking account of the reddening 
correction and the Ha luminosity density radiated from 



star-forming galaxies (Gallego et al. 1995), the fraction 
of AGN contribution to the total Ha luminosity density 
is about 15 % in the local universe. If we assume that 
the 15 % of the Ha luminosity density is radiated from 
AGNs, the corrected SFRD is 1.8±g;4 X 10~ 2 M & yr" 1 
Mpc" 3 . 

We note here that the error to psfr (and £(Ha)) is 
probably underestimated, since it does not include the 
effect of different correction methods and selection bi- 
ases. For example, adopting the different relation for 
correcting A^a may cause a different value of SFRD. 

We compare our result with the previous investigations 
compiled by Hopkins (2004) in Figure [H] We also show 
the evolution of SFRD derived from the observation of 
GALEX with mean attenuation of A™y s — 1.8, evalu- 
ated from the FUV slope (3 (fx oc A' 3 ) and the relation 
of ^4fuv = 4.43 -I- 1.99/3. If we adopt the more represen- 
tative value A™y = 1 (Schiminovich et al. 2005) deter- 
mined by using the i*d ust /Fuv ratio (Buat et al. 2005), 
their SFRD becomes smaller by a factor of 2, being sim- 
ilar to our SFRD. 

The left panel of Figure [6] shows the evolution of the 
SFRD as a function of redshift from z = to z = 2. 
The right panel of Figure [5] shows that as a function of 
the look-back time. It clearly shows that SFRD mono- 
tonically decreasing for 10 Gyr with increasing cosmic 
time. We note that the error of SFRD of our evaluation 
includes only random error, since we adopt the same as- 
sumptions as those in Hopkins (2004). 

Our SFRD evaluated above seems roughly consis- 
tent with but slightly smaller than the previous eval- 
uations, e.g., Tresse & Maddox (1998) and Fujita et 
al. (2003). Since we select emission-line galaxies with 
EW(Ha + [Nii]) bs > 12 A, our sample is considered to 
be biased against star-forming galaxies with small spe- 
cific SFR which is defined as the ratio of SFR to stellar 
mass of galaxy. Since our criterion is similar to that of 
Tresse & Maddox (1998) and the same as that of Fujita 
et al. (2003), we consider that the difference between our 
survey and theirs is not caused by the different criteria of 
EW(Ha + [Nn]) bs- As we mentioned in section 4.1, the 
SFRD of Fujita et al. (2003) was overestimated because 
of the contamination of [Oiii] emitters. On the other 
hand, the difference between Tresse & Maddox (19989 
and our work seems to be real; e.g., the cosmic variance. 

We discuss further the effect of the selection criterion of 
£W(Ha -I- [Nn]) obs > 12 A on the evaluation of SFRD. 
Being different from the previous Ha emission-line galaxy 
surveys using the objective-prism, the fraction of galaxies 
having EW(Ha) > 50 A is 12 % in our sample, which is 
similar to or less than the value of the local universe (15- 
20 %: Heckman 1998) and SINGG SRI (14.5 %: Han- 
ish et al. 2006). On the other hand, the fractions of 
the galaxies with EW(Ha) > 50 A are 42 % and 35 % 
in the KPNO International Spectroscopic Survey (KISS) 
(Gronwall et al. 2004) and UCM objective-prism sur- 
veys, respectively. Our sample seems to be not strongly 
biased toward galaxies with high equivalent width. Han- 
ish et al. (2006) showed that 4.5 % of the Ha luminosity 
density comes from galaxies with EW(Ha) < 10 A in 
local universe. If the fraction (4.5 %) is valid for the 
star-forming galaxies at z k 0.24, our estimate of SFRD 
would be about 5 % smaller than the true SFRD. 
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4.3. Spatial Distribution and Angular Two-Point 
Correlation Function 

Figure [7] shows the spatial distribution of our 980 Ha 
emitter candidates. There are some clustering regions 
over the field. To discuss the clustering properties more 
quantitatively, we derive the angular two-point correla- 
tion function (ACF), w(0), using the estimator defined 
by Landy & Szalay (1993), 

wm _ oow -m»w + m w (16) 

where DD(0), DR(9), and RR(6) are normalized 
numbers of galaxy-galaxy, galaxy-random, random- 
random pairs, respectively. The random sample con- 
sists of 100,000 sources with the same geometrical con- 
straints as the galaxy sample. Figure 0] demonstrates 
that our Ha emitter sample is quite incomplete for 
log£(Ha)(erg s^ 1 ) < 39.8. We therefore show the ACF 
for 693 Ha emitter candidates with log£(Ha)(erg s _1 ) > 
39.8 in Figure M The ACF is fit well by power law, 
w(0) = 0.013tEJ ™?e»-°- 88±0 03 . Recently, the departure 
from a power-law of the correlation function is reported 
(Zehavi et al. 2004; Ouchi et al. 2005). Such departure 
may be interpreted as the transition from a large-scale 
regime, where the pair of galaxies reside in separate ha- 
los, to a small-scale regime, where the pair of galaxies 
reside within the same halo. We find no evidence for 
such departure in our result. We however consider that 
the number of our sample is too small to discuss this 
problem. 

For Lyman break galaxies, brighter galaxies (with 
a larger star formation rate) tend to show more 
clustered structures than faint ones (with a smaller 
star formation rate) (e.g., Ouchi et al. 2004; 
Kashikawa et al. 2006). We also show the 
ACF of Ha emitters with larger Ha luminosity 
[logL(Ha)(ergs- 1 ) > 40.94 = log(0.1L*)] and that 
with lower Ha luminosity (39.8 < log L (Ha) (erg s _1 ) < 
40.94) in Figure M Both ACFs are well fit with 
a power law form: w(d) = 0.019±g;^gt6'- 1O8±0 - 05 
for objects with log L(Ha)(erg s _1 ) > 40.94, while 
w{9) = 0.01lt°;^6»-°- 84±0 - 05 for objects with 39.8 < 
logL(Ha)(erg s~ x ) < 40.94, respectively. We conclude 
that galaxies with a higher star formation rate are more 
strongly clustered than ones with a lower star formation 
rate. This fact is interpreted as that galaxies with a 
higher star formation rate reside in more massive dark 
matter halos, which are more clustered in the hierarchical 
structure formation scenario. 

It is useful to evaluate the correlation length ro of the 
two-point correlation function £(r) = (r/ro)~ 7 . A cor- 
relation length is derived from the ACF through Lim- 
ber's equation (e.g., Peebles 1980). Assuming that the 
redshift distribution of Ha emitters is a top hat shape 
of z = 0.242 ± 0.009, we obtain the correlation length 
of ro = 1.9 Mpc. Therefore, the two-point correla- 
tion function for all Ha emitters is written as £(r) = 
(r/1.9Mpc) -188 . The correlation length of Ha emitters 
with log L (Ha) (erg s^ 1 ) > 40.94 is 2.9 Mpc, while that 



of Ha emitters with 39.8 < logL(Ha)(erg s _1 ) < 40.94 
is 1.6 Mpc. These values are smaller than those evalu- 
ated for nearby L* galaxies (~ 7 Mpc) (Norberg et al. 
2001; Zehavi et al. 2005) and z ~ 1 galaxies (~ 4 - 5 
Mpc) (Coil et al. 2004). 

It is known that the correlation length is smaller for 
fainter galaxies in the nearby (Norberg et al. 2001; Ze- 
havi et al. 2005) and the z ~ 1 universe (Coil et al. 
2006). Figure [9] shows the relation between the L(Ha) 
and i?c _ band absolute magnitude Mr for our sample. 
Our sample includes many faint (Mr > —18) galax- 
ies. However, the correlation length for galaxies with 
-18 < M r < -17 (3.8 Mpc: Zehavi et al. 2005) is still 
larger than that of our sample. This discrepancy may 
imply a weak clustering of emission-line galaxies. 

5. SUMMARY 

We have performed the Ha emitter survey in the HST 
COSMOS 2 square degree field using the COSMOS of- 
ficial photometric catalog. Our results and conclusions 
are summarized as follows. 

1. We found 980 Ha emission-line galaxy candi- 
dates using the narrow-band imaging method. The 
Ha luminosity function is well fit by Schechter func- 
tion with a = -1. 35±o i3, logfc = -2.65tg;|£, and 
logI/*(erg s _1 ) = 41.94+0 1 ; 3 , 8 . Using the parameter set 
of Schechter function, the Ha luminosity density is eval- 
uated as 2.7j:°;g x 10 39 erg s _1 Mpc" 3 . If we adopt 
the AGN contribution to the Ha luminosity density 
is 15 %, we obtain the star formation rate density of 
1.8+q;4 x lO" 2 M yr- 1 Mpc- 3 . This error includes only 
random error. Our result supports the strong increase in 
the SFRD from z = to z = 1. 

2. We studied the clustering properties of Ha emitters 
at z ~ 0.24. The two-point correlation function is well fit 
by power law, w(9) = 0M3± o Q a o o o l9-° S8±om , which leads 
to the correlation function of (r/1.9Mpc)~ 188 . We can- 
not find the departure from a power law, which is recently 
found in both low- and high-z galaxies. Although the 
power of —1.88 is consistent with the power for nearby 
galaxies, the derived correlation length of ro = 1.9 Mpc 
is smaller than that for nearby galaxies with the same 
optical luminosity range. This discrepancy may imply a 
weak clustering of emission-line galaxies. The galaxies 
with higher SFR are more strongly clustered than those 
with lower SFR. Taking account of the fact that the SFR 
of a luminous galaxy is higher than that of a faint galaxy, 
this result is consistent with the fact already known that 
the luminous galaxies are more strongly clustering. 
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We would also like to thank the Subaru Telescope staff 
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Fig. 1. — Diagram of iz — NB816 vs. ./V.B816 for all objects classified as galaxies in the ACS catalog. The horizontal solid line corresponds 
to iz — NB816 = 0.1. The dashed lines show the distribution of 3<r error, the dot-dashed line shows the limiting magnitude of iz. Since the 
total i'-magnitudes of galaxies in the official photometric redshift catalog are brighter than 25, iz magnitudes of most of them are brighter 
than the limiting magnitude. 
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Fig. 2. — Diagrams between B — V vs. V — r'. Top: Colors of model galaxies (CWW) from z = to z = 3 are shown with dotted 
lines: red, orange, green, and blue lines show the loci of E, Sbc, Scd, and Irr galaxies, respectively. Colors of z = 0.24, 0.64, 0.68, and 1.18 
(for Ha, [O in], H/3, and [O II] emitters, respectively) are shown with red, green, light blue, and blue lines, respectively. Orange asterisks 
show Gunn and Stryker (1983)'s star. Bottom: Plot of B — V vs. V — r' for the 6176 sources found with emitter selection criteria. In this 
diagram, Ha emitters are located above the black line, that is adopted by us as one of the criteria for the selection of Ha emitters. The 
980 Ha emitters are shown as black dots and other emission-line galaxy candidates are shown by gray dots. Galaxies in GOODS-N (Cowie 
et al. 2004) with redshifts corresponding to Ha emitters, [Oni] emitters, H/3 emitters and [On] emitters are shown as red, green, light blue 
and blue open squares, respectively. 
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Fig. 3. — Diagrams between B — r' vs. i' — z' . Top: Colors of model galaxies (CWW) from z = to z = 3 are shown with dotted 
lines: red, orange, green, and blue lines show the loci of E, Sbc, Scd, and Irr galaxies, respectively. Colors of z = 0.24, 0.64, 0.68, and 1.18 
(for Ha, [O III], H/3, and [O II] emitters, respectively) are shown with red, green, light blue, and blue lines, respectively. Orange asterisks 
show Gunn and Stryker (1983)'s star. Bottom: Plot of B — r' vs. %' — z' for the 6176 sources found with emitter selection criteria (black 
dots). In this diagram, Ha emitters are located above the both of black lines, that is adopted by us as one of the criteria for the selection 
of Ha emitters. The 980 Ha emitters are shown as black dots and other emission-line galaxy candidates are shown by gray dots. Galaxies 
in GOODS-N (Cowie et al. 2004) with redshifts corresponding to Ha emitters, [Om] emitters, H/3 emitters and [On] emitters are shown as 
red, green, light blue and blue open squares, respectively. 
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Fig. 4. — Our Ha luminosity function (filled squares and thick solid line) and Ha luminosity functions in previous works. The Tresse 
& Maddox (f998)'s Ha luminosity function at z < 0.3 is shown with the dashed line. The Ha luminosity functions derived by Fujita et al. 
(2003), Hippelein et al. (2003), and Ly et al. (2007) are shown with the dotted line, the dot-dashed line, and dashed and double-dotted 
line, respectively. Data points of Ly et al. (2007) 's Ha LF are shown as gray crosses. 




Fig. 5. — B — Rc vs. Rc — Ic color - color diagram of model galaxies. Colors of z = 0.24, 0.64, 0.68, and 1.18 (for Ha, [O III], H/3, 
and [O II] emitters, respectively) are shown with red, green, light blue, and blue lines, respectively. The loci calculated by using GISSEL96 
(Bruzual & Chariot 1993) are shown by solid lines and those calculated by using CWW are shown by dashed lines. Galaxies in GOODS-N 
(Cowie et al. 2004) with redshifts corresponding to Ha emitters, [Om] emitters, H/3 emitters and [On] emitters are shown as red, green, 
light blue and blue open squares, respectively. Fujita et al. (2003) selected galaxies above the black solid line as Ha emitter candidates. If 
we reselect Ha emitter candidates as sources above the black dashed line, some of the Ha emitter candidates (black dots) do not satisfy 
the new criterion. 
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Fig. 6. — Star formation rate density (SFRD) at 2 rj 0.24 derived from this study (large red filled circle) shown together with the 
previous investigations compiled by Hopkins (2004). SFRDs estimated from Ha, [Oil], and UV continuum are shown as orange open circles 
(Perez-Gonzalez et al. 2003; Tresse et al. 2002; Moorwood et al. 2000; Hopkins et al. 2000; Glazebrook et al. 1999; Yan et al. 1999; Tresse 
& Maddox 1998; Gallego et al. 1995), green open diamonds (Teplitz et al. 2003; Gallego et al. 2002; Hogg et al. 1998; Hammer et al. 
1997), and blue squares (Wilson et al. 2002; Massarotti et al. 2001; Sullivan et al. 2000; Cowie et al. 1999; Treyer et al. 1998; Connolly et 
al. 1997; Lilly et al. 1996). The light blue open squares show SFRDs based on the UV luminosity density by Schiminovich et al. (2005), 
assuming Apyy = 1.8. An orange open square and an orange open diamond show SFRDs at z ss 0.24 derived by Fujita et al. (2003) and 
Ly et al. (2007), respectively. In the left panel, we show the evolution of SFRD as a function of redshift, and in the right panel, we show 
it as a function of lookback time. 




— 1 1 1 1 1 — 1 — 1 — 1 — i 1 1 1 1 1 1 — I I 11 

150.5 150 149,5 150.5 150 149.5 

RA (degree) RA (degree) 

Fig. 7. — Spatial distributions of our H a emitter candidates (black filled circles and black dots). Gray open squares in the both panels 
show our survey area. The shadowed regions in the right panel show the areas masked out for the detection. We show the luminous Ha 
emitters [log L(Ha)(erg s -1 ) > 40.94] as large filled circles and the faint Ha emitters [39.8 < logL(Ha)(erg s -1 ) < 40.94] as small filled 
circles. Ha emitters with logL(Ha)(erg s _1 ) < 39.8 are shown as black dots. 
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Fig. 8. — Angular two-point correlation function of all H a emitter candidates (filled circles), bright H« emitter candidates 
(log L(Ha)(crg s _1 ) > 40.94: open squares), and faint Ha emitter candidates (39.8 < log L(Ho)(crg s _1 ) < 40.94: open triangles). 
Solid line shows the relation of w(9) = O.O130 -0 88 . Dashed line shows the best-fitting power law for bright ones, w(6) = O.O190 -1 - 08 , and 
dotted line shows that for faint ones, w(9) = 0.0116 -0 ' 84 . 




